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Abstract 

We discuss general implications of the local spin-triplet pairing among 
fermions induced by local ferromagnetic exchange, example of which is the 
Hund's rule coupling. The quasiparticle energy and their wave function are 
determined for the three principal phases with the gap, which is momentum 
independent. We utilize the Bogolyubov-Nambu-De Gennes approach, which 
in the case of triplet pairing in the two-band case leads to the four-components 
wave function. Both gapless modes and those with an isotropic gap appear 
in the quasiparticle spectrum. A striking analogy with the Dirac equation is 
briefly explored. This type of pairing is relevant to relativistic fermions as 
well, since it reflects the fundamental discrete symmetry-particle interchange. 
A comparison with the local interband spin-singlet pairing is also made. 
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I. Introduction 



The discovery of superconductivity in the orbitally degenerate system Sr2Ru04, which 
is closely related to both ferromagnetic SrRuO^ and antiferromagnetic and Mott insu- 
lating [§] Ca2Ru04, poses a question about the role of short-range Coulomb and exchange 
interactions in stabilizing the spin-triplet superfluid state [Q]. In the case of Mott-Hubbard 
insulators the kinetic exchange interaction 0] plays an essential role in stabilizing antiferro- 
magnetism. This interaction is also instrumental in the form of real space pairing |^ in driv- 
ing the system close to the Mott-Hubbard boundary towards spin-singlet superconducting 
state. In the case of orbitally degenerate systems, the ferromagnetic and antiferromag- 
netic kinetic exchange interactions compete with each other for the number of electrons 
per atom n > 1. Ferro magnetism (with a possible orbital ordering) usually wins for 
n —>■ 1, whereas the antiferro magnetism takes over when n ^ d, where d is the orbital 
degeneracy. This type of competition should also be present in Sr2Ru04, in which 4d^ con- 
figuration of Ru*~^ contains two holes in t2g shell composed of nominally triply degenerate 
de = {dxy,dyz,dzx) orbitals. The two-dimensional antiferromagnetic spin fluctuations have 
been indeed observed in Sr2Ru04 system . From the symmetry point of view d^y does not 



mix with dyz and d^x, so the fluctuations can be ascribed |I0|] as solely due to the electrons 
in dxy band. The Hund's rule coupling between dxy and the remaining two bands {dyz, dzx) 
must than suppress the formation of the antiferromagnetic state. In effect, we are left with 
two electronic liquids: the doubly degenerate and hybridized dyz — dzx band containing ap- 
proximately one hole and the d^y band containing the other. It must be underlined that 
all t2g holes are delocalized, since one observes a well defined Baber-Landau-Pomeranchuk 



T ) contribution to the resistivity in both x — y {RUO2) plane and in c direction [p 



From what has been said above it is important to formulate first the model of local 
pairing represented a doubly degenerate (or almost degenerate) band coupled by the Hund's 
rule and characterize the possible spin-triplet solutions induced by the Hund's rule (ferro- 
magnetic) exchange. This type of model has been formulated by us recently ||12|. We have 



shown there that sizeable (of the order of bare bandwidth) Coulomb correlations renormal- 
ize the system properties, i.e. lead to an almost localized Fermi liquid with a nonretarded 
real-space and spin-triplet pairing. The renormalized Fermi-liquid nature of our fermionic 
system will be a starting point in this paper, in which we consider basic features of the 
superconducting state such as the quasiparticle wave function (in the Fock space) and their 
energies. We list the possible solutions for our effective model with interorbital pairing. 
The question of coexistence of the Al state with ferromagnetism, as well as the competition 
with the orbitally ordered-spin ferromagnetic state has been discussed separately [|^. We 
believe that the present two-band model stands on its own ground, independently of the 
detailed nature of Sr2Ru04 superconductivity (which should include the third band and 
the anisotropic interband hybridization) and must be considered separately, to amplify the 
physical plausibility of this mechanism of spin-triplet pairing (see also the discussion at the 
end). This is particularly so because the Hund's rule and associated with it ferromagnetic 
fluctuations |jT^ represent probably the most natural determinants of spin-triplet pairing 
under these circumstances. Also, the present real-space pairing |13] represents is formally 
analogous to the spin-singlet pairing |^ and additionally, reflects a fundamental symmetry 
- the particle interchange. So, it contains fundamental physics in the sense, that the nature 
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of the ground state, i.e. that of the spin-triplet superconductor, can appear instead of or 
together with an itinerant ferromagnetism. 



II. Nambu-De Gennes method for the triplet pairing in 
the two-band case 

We consider a degenerate two-band Fermi-hquid system coupled by a local triplet pairing. 
The corresponding effective Hamiltonian is of the simple form ||T^ 

= ^ Ekiakia^kia - 2 J ^ Aj^Aim, (1) 

kcr/=l,2 im 

where i?kz are the quasiparticle energies with enhanced masses by the band narrowing factor 
(calculated self-consistently in the bands / = 1,2, J ~ Jt^ is the effective Hund's 
rule coupling (the local interorbital exchange), and t"^ is the probability of having interorbital 
local spin-triplet configurations, characterized by the creation operators A[ = al^ali,p = 

^lil^li'l^ and Aq = -^{a\i^ali,^ + a\i^ali,-^) for / 7^ The local exchange origin of the second 
term derives from the exact relation between the pairing operators in real space and the 
full exchange operator projecting the corresponding two-particle state onto the spin-triplet 
configuration. 

1 3 
^ Al^Aim = Sil ■ Sii> + -nuriiv, (2) 

m=-l ^ 

where Sj/ and rin are respectively the spin and the particle number operators for electron on 
site i and orbital /. Explicitly nu = J2a ^ua, nua = alig-aua, whereas the spin operators Su = 
{Su^S~i,S^i) = (a|^|a.^,aj^a.^|, (l/2)(n.^| - n.^)). The right-hand side of (D represents 

thus the full exchange operator. 



After making the BCS-type approximation in the local form ||I4 

-^im^im — -^im ^ -^im ^ ~l~ ^ -^im ^ -^im ^ -^im ^ (3) 

we can cast Hamiltonian (|I]) into the four- component form, which in the reciprocal (k) space 
takes the form ||12 

HbC5 = EfkHkfk + E^k2, (4) 

k k 

where the corresponding Nambu operators take the form: f,^ = (/j[^|) /kij, /_k2|' ^-k2 
fk = (fk)^ and the Hamiltonian matrix for selected k state reads 



Hi 



/ ^kl - /i, 


0, 


Ai, 


Ao \ 


0, 


Eki — /i. 


Ao, 


A-i 


A*, 


A* 


— -Ek2 + 





V AS, 


A* 


0, 


-Ek2 + /i / 



EkiOTQ, A 
A*, -Ek2ao 



(5) 
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where (Tq = 1 is the unit 2x2 matrix, and /i is the chemical potential. The superconducting 
gap is parametrized as = — 2JI]k < fliaf-k2a' with m = ((j + cr')/2, and a, a' = ±1. 



The 2x2 matrix A is parametrized in the usual form [16 



A = i(d ■ a)(T,, 



-drc + id. 



yi 



dz 

dj^ I tdf. 



(6) 



where a is composed of the three Pauli matrices, whereas the vector d in spin space has 
the components dx = (A_i — Ai)/2, dy = (A_.i + Ai)/2, and dz = Aq. The form 
is a generalization of the Nambu representation to the triplet case with three, in general 
different, gaps A^- 

It is straightforward to introduce the 4x4 Dirac matrices 



/3 



1, 
0, 







and a,- 



0, (Xi 
0",;, 



and then rewrite (|) for the degenerate case i?ki = -Ek2 and for A^ = Aj^ in the form 

Hk = /9(^k-/i)+^(d-5)S2, (7) 

where 

""'-W J 

is the y component of the relativistic spin operator. We discuss in detail the simple situation 
of degenerate electrons (i?ki = -E'k2) with a real gap A^ in the next Section. 

One can also look at the approach from a different prospective. Let us introduce the four 
component wave function for a single quasiparticle in the suprconducting phase propagating 
in the real space as follows 



^(x,t) 



E 



/ ^ik/kit 



k2t 
^4k/!,,| / 



exp 



n 



(8) 



where T/'^k are the quasiparticle amplitudes which are determined for each eigenstate (see 
below). In this representation the Bogolyubov-De Gennes equation for a single quasiparticle 
in the superconducting states reads: 



thdt^ = /3{^k(k -) 



^}^ + z(d-a)S2^, 



(9) 



where £^k(k =^ represents now the differential operator (l/z)V replacing the wave vector 
k in the dispersion relation Ey_ for quasiparticles. In the effective-mass approximation and 
in the stationary case this wave equation for quasiparticles in the superconducting phase 
has the following form 



4 



A 



^3 



2m* 



( ^1 \ 

^2 



^ Ai^3 + AoV'4 \ 

Ao^As + A_ii/'4 
AiT/^i + Ao^2 

V Ao^/'i + A_i^2 y 



(10) 



where = 'i/'^(x) and A is an eigenvalue of quasiparticle state in the superconducting state 
with the above 4-component wave function (A^ are regarded as real). The validity of this 
equation goes beyond the simple solution (|), as one can include the magnetic and electric 
fields and other inhomogeneities if they appear on the mesoscopic or macroscopic scale. In 
the next Section we will use explicitly the momentum representation of Eqs.(|r^), as we 
will discuss exclusively homogeneous superconducting states. We will return to Eqs.(|10|) 
when discussing the general features of this Hamiltonian in Section IV. One should also 
note that finding the eigenvalues for Hamiltonian in the forms or (J^) can be achieved by 
diagonalizing of the matrix 4 x 4 in general discussed in analytic terms in Appendix 

A. 



III. Superconducting states and their quasiparticles 

We now discuss three principal solutions of Eq. (1^) by taking ^/'^(x) = ip^ exp(ik-x) / \/V, 
where V is the system volume. We also assume that A^ = A*, (e.g. neglect the applied 
magnetic fields), since we consider only spatially homogeneous solutions. Namely, rewriting 
Eq.(|TUp in components we obtain the combinations 



and 



Ki^i + ^2) = (^k - /i)(V^i + ^2) + (Ai + Ao)^3 + (Ao + A_i)V^4 

A(^3 + ^4) = -(^k - yU)(V'3 + ^4) + (Ai + Ao)^i + (Ao + A_i)^2, 



A(^/'i - V2) = (^k - yu)(^i - ^2) + (Ai - Ao)^3 + (Ao - A_i)V'4 
X{ij3 - i^i) = -(^k - /i)(^3 - ^4) + (Ai - Ao)^i + (Ao - A_i)^2. 



(11) 



(12) 



Such combinations of particle {ipi and tp2) and hole {ips and ipi) components contain basic 
symmetry, as we will see on example of particular solutions, which we discuss next. 

A. Isotropic solution: Aq = A 1 = Ai = A 

In that situation Eqs.(pJ]) - ([1^) take a simple form 

(13) 



f H^i + ^2) = (^k - fi'KA + ^2) + 2A{^3 + ^4) 
\ A(^3 + i^d = -(^k - /i)(^3 + ^4) + 2A(^i + ^2) 



and 



\ A(^3 - V'4) - 



(Ek - fl){iJl - tp2) 
-(Ek - /i)(V^3 - ^4)- 



(14) 



5 



The first two equations lead to the modes with a gap 



(15) 



For those two modes ipi = ip2 and ips = and their eigenstates are characterized by the 
following quasiparticle operators 



"k = Mk^ (/nt + /, 



and 

with the Bogolyubov coherence factors 



1 



'"■75 l^--T 



-k2 



k2l 



(16) 



(17) 



72 V 



^^k 



1 

71 



Al 



(18) 



The quasiparticle operators contain symmetric combinations {f^^^i^ + f^^-^^)/ V2 and (/_i(.2| + 
/ , I )/V2. The wave function is symmetric with respect to particle-spin interchange (t ^ 



I) and describes quasiparticle states of energy ±Ak, respectively. 
Eqs.(|Tip lead to the gapless modes of the form 

A = Ak3,4 = ±(-E'k — /i), 
and correspond to the eigenstates characterized by the operators 



7k 



and 6^ 



y/2 V'^-k2t •^-k2i 



(19) 



(20) 



and constitute the antisymmetric-in-spin operators, representing the unpaired electrons. 
These gapless modes disappear when the gap components are not equal, as shown in Ap- 
pendix A. One should note that the gapless modes appear even though the superconducting 
gap here is k-independent. 

Combining the solutions (0) - (p!8| ) and (|19| - ( pO]) we can express the original ("old") 
particle operators in terms of quasiparticle ("new") operators in the following manner 



/ /kit ^ 

•'-k2t 
V ^-k2i / 





f 


Mk, 


Vk, 


1, 




/ 


ttk 


\ 


1 


Mk, 


Vk, 


-1, 







/51k 










Uk, 


0, 


1 




7k 






[ 






0, 


-1/ 


V 


5lk 


/ 



(21) 



This transformation is necessary for determining the self-consistent equation for the gap and 
for the chemical potential /i. First, we rewrite the Hamiltonian (11) in the diagonal form 



H 



E ^k(aL«k - /3-k/3lk) + ^k(7l7k - 5-k5lk) + E ^k 

k k 



6 



E Ak(4«k + /5lk/5-k - 1) + £^k(7i7k + 5^k5-k). (22) 



The equation for the gap e.g. Ai =< /j[^|/^j^2| ^ obtained by substituting the relevant 
transformed operators in (^) to Ai. In effect, we obtain the usual BCS form = E\^ — ^) 



+ + 1 A //3Je2+4A2\ 

where /? = {kBT)^^ . So, the gap equation has two solutions: 1" A = 0, 
2° 



The last equation tells us that the physical gap is 2 A. The self-consistent equation for the 
chemical potential must include gapless modes, i.e. takes the form 

1 2 
= T7 E < fliahia + /k2CT/k2a >= T7 E < "k"k + 7k7k > • (25) 

k(T k 

Normally, as we shall see, A <^ and hence approximately half of all particles will have 
the spectrum gapped. The details must be analysed numerically for a concrete structure of 
the density of states. In the limit ly -C J we have the estimate of the gap value at T = 
in the form A = iW/2) exp(-W'/(2 J); this yields the value A/W" ~ lO'^ - IQ-"^, or in the 
regime 1 — lOK for W ~ leV and J ~ O.ll^. 

We need also the expression for the ground state energy, as various solutions are possible. 
In the present case, this energy can be written as 

^ = 1 Ei^El + A^ < aia^ > < 7^k > -^JeI + A^} + ^, (26) 
where A = 2A . 



B. Equal-spin pairing: Aq = 

To obtain the explicit solution we now combine separately the first and third components 
of Eq. ( p!OD on one side, and the second and the fourth on the other. Adding and subtracting 
the corresponding terms we obtain: 

r (Ai - A)(V^i + V^s) + (^k - - ^3) = . . 

\(Ek-/i)(V^i + ^3)-(Ai + A)(V^i-V^3) = 0, ^ 

and 

r (A_i - A)(7/>2 + V^4) + (E^ - /i)(^2 - ^4) = . . 

1 (Ek - fi){^2 + ^a) - (A_i + X){^2 - ^4) = 0. ^^""^ 



Thus, the two pairs of components (^) and ( pSD separate from each other and it is sufficient 
to solve e.g. the first system (^71) to be able to reproduce the other. Explicitly, the two 
solutions can be combined into the form, in which the eigenvalues take the form 



A = Aki...4 = ±V(^k - + = ±A 



(29) 



where for each spin orientation cr = ±1 of the quasiparticles we have two solutions with 

the gap iy' (-Ek — yu)^ + ^a- The quasiparticle operators {a^a, P-ka) diagonalizing Hamil- 
tonian (Hf) in this case are: 



M^^;^ (/kia + /. 



t 



■k2cr 



i^) ^ (f 



/I 



and 



-kcr 



with the coherence factors 



klO" 



+ /1 



k2cr 



1 



A, 



1/2 



A 



1 

71 



klO" 



k2cr ; ) 



/-k2(T ; ) 



A 



1/2 



A 



(30) 



(31) 



(32) 



In general, we have two gaps Aq- = (Ai, A_i). In the situation Act = A_cr = A we have a 
doubly (spin) degenerate solutions. It can be easily verified that the operators (^) and (^Tj) 
obey the fermion anticommutation relations. The diagonalized Hamiltonian has the form 



^ = E («k(T"ka + 4a/5ka - l) + ^0 



(33) 



To determine the gap equation we have to find the transformation which is reverse of 
(BO) and (P). It is of the form 



/ /kit \ 

4k2t 

V^~k2iy 



1 

71 



u 



(+) 



0, 



(+) .,(+) 

k ' 



..{+) 



u 



-u 



(+) 



0, 

k + '^k ' 
0, 







u 



(-) 



,(-) 



V 



0, 



0, 



u 



{-) 



(-) 

k ) 



(") 



«kt 

^IkT 
"ki 



(34) 



The difference with the isotropic pairing (^) is that here the coherence factors appear in 
combinations. Those appear also in the self-consistent equation for the gap 



< /-k20"/kia 



>= -iu 



(0-)2 



)[1- < Oko-ako- > - < /3ko-/5ko- >] 



(35) 



In result, the self-consistent equation will have the following three solutions: 1° A^r = 0,2° 
A(_(j) = 0, but Act 7^ is the solution of equation 



AT ^ 



k v/(^k - /i)2 + A2, 



tanh^J(Ek-/i)2 + A2, 



(36) 
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3° Act 7^ 0, A_o- 7^ 0, and each of them is determined from equation (|36|). 

One should note that the couphng constant above (J) is the same as for the isotropic 
phase (cf. Eq.(p^). For the sake of completness, we reproduce the ground-state-energy 
expression which is 

1 



- ^ E ^(Ea - /i)^ + + - E E^. (37) 

This phase represents a starting point when discussing the coexistence of ferromagmnetism 
and the spin-triplet superconductivity. 



C. Spin-polarized phase: Aq = Aj = 

In this limit the system is totally spin polarized, i.e. is a spin- saturated supercon- 
ductor. In that limit we recover again the spectrum both with and without gap, i.e. 

Aki,2 = ±y''(-Ek — /^)^ + A|, Ak3,4 = ±(-E'k — A^)- Thus paired and unpaired states coex- 
ist also in this phase, as can be easily seen from Eqs. (|30D - (|3lD , which yield the form 
written there for cr = t and Oj^j = f.^^^ and /^^,^| = ~fl^2[' 

Summarizing cases A-C, the lowest energy will have the homogeneous state with A| = 
A| = Aq so that the effective gap is equal to 2A. The most interesting feature of the results 
is that the gapless modes coexist in cases A and C and represent half of the spectrum. 
Also, the condensed phases described by the cases A-C above correspond roughly to the 
solutions for superfluid ^iJe, which are labelled B, A, and Al. However, under the present 
circumstances here we have momentum independent gaps, since the pairing is of the local 
(intrasite, but interorbital) nature. For the sake of comparison we present in Appendix B 
the case of spin-singlet pairing induced by the same type of local interband pairing induced 
by antiferromagnetic exchange. 



IV. Remark on the triplet pairing for relativistic 
fermions 



The two-band situation with a local ferromagnetic exchange can be easily generalized 
to the explicitly relativistic form modelling thus the triplet configuration of spin, isospin or 
color (the singlet case was considered by Nambu and Jona-Lasinio |]l6l and in [jl9l)- The 
paired quasiparticles PD| obey the following modified Dirac wave equation 



ihdt^ = {ca-p + Pmc^)'^ + i{d ■ a)S2^, 



(38) 
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where the last term supplements the Dirac equation with the contact pairing. By taking 



the analogy with the original approach by Nambu |^ one can write down the stationary 
version of this equation as the following system in the two-component (Weyl) representation 




{ca ■ p 




+ mc 




Ao^3 + A_iV^4 



(39) 



A 



^3 
^4 



[ca 




+ mc 




(40) 



This system of equations can be directly compared with Eq.(|TO]) for nonrelativistic electrons. 
In the present situation the mass term mixes the upper and the lower two components of 
the bispinor, as does the pairing part. Apart from a modification in the kinetic-energy 
part, the system of equations (^)-(|40|) can be solved in the manner as discussed in Section 
III. The detailed discussion must include also the gauge fields, which lead to one-exchange 
pairing (analogous to the phonon-mediated pairing), a topic intensively discussed in recent 
literature HI. In general, the singlet pairing iUHl, |19[ - ||2^ is mutually exclusive with the 



^ . In general, the singlet pairing [|n 

triplet pairing proposed here and therefore, the latter requires first a detailed discussion of 
exchange interactions represented by relativistic spin {Sj}. 



V. Discussion 

In this paper we have formulated the quasiparticle language for local triplet pairing 
between fermions (interband pairing in the nonrelativistic case) induced by the local (Hund's 
rule or Dirac) exchange. In particular, we have determined explicitly the quasiparticle states 
and the De Gennes wave equation for them, which can be useful when considering spatially 
inhomogeneous situations |2^. The principal feature of the results is the existence of the 



gapless modes, existence of which can also be proved on a phenomenological level |^ . The 



circumstance that the pairing is induced by the ferromagnetic exchange means that this 
interaction can lead not only to an itinerant ferromagnetic state, but also to either spin- 
triplet superconductor or to a coexistence of both these states (for a brief discussion of 
this issue see Ref. [0). The present paper represents only the first step in this direction. 
Furthermore, our mechanism of pairing expressing the fundamental symmetry (the particle 
interchange) may have an important astrophysical application: the pairing in the neutron- 
proton matter in pulsar, but this intriguing possibility requires a separate study. 

Two problems should be tackled next. First, the analysis of the Meissner effect, since 
in the present situation the orbital diamagnetism will compete with the ferromagnetic spin 
polarization (particularly, if the triplet superconducting and ferromagnetic phases can co- 
exist). An intriguing question here is: can we reach the limiting superconducting phase 
(corresponding to Al phase in the case of superfluid ^He), the critical temperature of 
which can be enhanced by the applied magnetic field? 

Second, one should derive microscopically the Ginzburg-Landau equation for the con- 
densed pairs. Note that the De Gennes equation (^ or ( pUf ) is useful in describing the 



quasiparticle tunneling, whereas the Ginzburg-Landau equation is useful when consider- 
ing the Josephon (pair) tunneling. Here an intriguing question to what extent the gapless 
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quasiparticles influence the tunneling between the spin-singlet and the spin-triplet super- 
conductors or between the triplet superconductor and normal metal. We should be able to 
see the progress in answering those questions in a near future. 

Finally, returning to the question of the nature of the paired state in Sr2RuOi one can 
make the following two remarks. First, the existence of gapless modes in B- and Al-like 
phases leads to the persistence of the linear term in the specific heat in the superconduncting 
phase at its 50% value, if the pairing is the pure spin-triplet state of electrons pairs derived 
from bands. The recent measurements [^] in very pure samples contradict such earlier 
claims |]TT| that a half of the linear specific heat survives in the superconducting phases. 



Does that mean that the full phase diagram involves more than one phase depending on 
the doping degree, as in the heavy-fermion system Ui-xTh^Beis ? In connection with 



this one can say that because of the reasons mentioned in Section I it is conceivable that 
a singlet pairing in d^y band induced by antiferromagnetic fluctuations can compete in 
the triplet state in the other two bands The nature of the resultant state should be 
determined then. We should be able to see a progress in those matters in near future. 

Second, an important question concerns the nature of the pairing potential. In more 



standard approach [^], |jT^ one introduces the effective triplet pairing via the paramag- 



non exchange. In that situation the coupling constant is determined by the susceptibility 
x(q = k — k') and hence, is wave vector (q) dependent. In the approach developed here 
the exchange interaction itself provides real space pairing, as in the case of high temper- 
ature superconductors 0. In the case of Hund's rule coupling the pairing potential is 
then k- independent. We believe that the latter approach is relevant when the particles are 
strongly correlated. Sr2Ru04^ is a systems close to (but below) the Mott-Hubbard local- 
ization threshold, i.e. the halfway between the weak- correlation and the strong-correlation 
asymptotic regimes. Therefore, the real space pairing is certainly worth of analysing, as it 
allows for an analytic approach. 

One methodological remark at the end. In the analytical analysis of the spin-triplet 



pairing one usually uses fl^, p9[ the d vector in expressing the pairing part. Here we 



decided to use the original BCS gap parameters, a completely equivalent procedure, but 
probably a bit more direct, at least in the spatially homogenous situation. In connection 
with this a difference of the present description with that for superfluid helium — 3 should 
be stressed. Namely, in the helium — 3 case, the L = 1 orbital moment 1 and the spin 
vector d determine the (many-component) nature of the gap. The d.l and dxl combinations 
determine the order-parameter dynamics. Here, no 1 vector appears and therefore, the order 
parameter can have up to 3 independent components. 
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Appendix A: General solution: Ai 7^ A2 7^ A3 7^ Ai, E'ki 7^ 



The most general case of finding the eigenvalues for quasiparticles in the superconducting 
phase is to diagonalize matrix 4x4, i.e. to solve the equation (we assume that £'ki,2 — = -E'k 
and that Ai = A*) 



det 



1 -Eki — A, 


0, 


Ai, 


Ao 


\ 


0, 


-E'ki — A, 


Ao, 


A_i 


Ai, 


Ao, 


—E\a — A, 







V Ao, 


A-i, 


0, 


—E\a — 


A/ 



= 0. 



(41) 



By a straightforward evaluation one obtains the eigenvalues 



1 1 r ~ r ^1/2 

Aki..4 = -(^ki - ^^k2) ± - {E^i + E^^f + 2A2 ± 2VA4 - 



where 



is the total gap, and 



A = (A? + A^ + 2A^)i/2, 



(42) 



(43) 



A2-AiA_i|, 



(44) 



is its anisotropy in fermion-pair spin space. In this general case all four modes are each 
with a different gap and the results reduce nicely to the eigenvalues discussed A-C 
in main text. In general, the superconducting coupling at the level of energies amounts 
to hybridizing the different fermion fields (/ = 1,2) and their spin (cr t, J,) states. The 
general form of the eigenstates can also be found in a straightforward manner, but will not 
be discussed here. 
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Appendix B: Local spin-singlet pairing in two-band 

case 



For the sake of comparison with the spin-triplet case we outline here the solution for 
the corresponding spin-singlet situation. In this case the Hamiltonian with the spin-singlet 
exchange coupling has the form in the real space 



kcri=l,2 ill' 



(45) 



Note that now the exchange operator in the present situation differs from (0) introduced in 
the triplet case. Introducing the corresponding local pairing operators in the singlet state 



— ^ \ Qj -^Qj t Qj \ Qi ^ 
V2 V »lT ill i2] 



(46) 



we can write down the second term in 



as -2JY.iB\B,. 



After taking the space Fourier transform, and including only (k, — k) pairs we obtain 

J_ 

kcr kk' 



H = E E^.ln^.l<y - ^ E {fl]fL2[ ' ^[^^2]) (/-k2i/kiT " /-k2T/kii) • (47) 



Making subsequently, as in Section II the BCS approximation, we obtain (the chemical 
potential is included in E]^i = E]^i — ji) 



kia 



flif'r)+H.C.}-^N. 



klt''-k2i •'kl|''-k2t 



2J 



where 



A 



-jy E < /lllt^-k2i > • 



(48) 



(49) 



Introducing, as before, the four dimensional vectors = (/j[^|, /j^^' /-k2|' •^-k2|) their 
conjugate as one column vectors, we can write down the Hamiltonian in the form of the 
following 4x4 matrix 



/ -E^kl, 


0, 


0, 


A 


\ 


0, 


-Ekl, 


-A, 







0, 


-A*, 


—Ek2, 







V A*, 


0, 


0, 


— -Ek2 





/ /kit ^ 

•^-k2t 
V ^-k2i J 



(50) 



with Eq = 2 J2k E]fi2 + A^/ (2J). Diagonalization of this 4x4 matrix leads to the eigenvalues 

(51) 



A 



1,2 



2 (-^kl — Ek2) ± 



-E'ki + E] 



k2 



+ |A|2. 
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Both eigen modes are doubly degenerate and with an isotropic gap A. We take the form of 
usual dispersion relation for degenerate bands i?ki = -£'k2- 

The corresponding combinations of the wave function components are (for A = A*) 



and 



{Eki - A)(^i + ^2) - A(^3 - ^4) = 
A(V'i + iJ2) + (^k2 + A) (^3 - ^4) = 0, 



(Eki - A)(V^i - V^s) + A(V^3 + 7/'4) = 

A(^l - ^2) - (^k2 + A)(V^3 + V^4) = 0. 



(52) 



(53) 



For the sake of simplicity we consider here only the case Ey_i = E\^2 = Ey_, a.s it provides the 
main character of the eigenstates. Moreover, it is sufficient to consider only the system (|52D 
due to the double degeneracy of the eigenvalues. By standard method (including the wave 
function normalization) we obtain the quasiparticle operators 



"k = ^^k 



^ (/kit + /kii) + (/!k2T " ^-k2i 



V2 



and 



/3lk 



V2 



k2 



k2| 



(54) 



(55) 



where 



Mk 





Eu 



1/2 



El + A\ 



(56) 



Again, we have a combination of the two types of states. The corresponding wave equation 
which replaces the Bogolyubov-De Gennes equation for one-band singlet superconductor 
reads in the effective-mass approximation 



ihdt 



^3 



2m' 



-V^ -/i 



( ^ 




( ^ 


^2 


+ A 


-V^3 










\ i^i J 



(57) 



We see that the pairing couples explicitly the particle and hole components ( ipi with ^4, 
iIj2 with —ip3, etc.). This equation forms a basis for the discussion of inhomogeneous paired 
states. 
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